arXiv:1503.05749vl [astro-ph.SR] 19 Mar 2015 


Astronomy & Astrophysics manuscript no. arXbcep ©ESO 2015 

March 20, 2015 


On X-ray pulsations in p Cephei-type variables * 

(Research Note) 

L. M. Oskinova 1 , H. Todt 1 , D. P. Huenemoerder 2 , S. Hubrig 3 , R. Ignace 4 , W.-R. Hamann 1 , and L. Balona 5 

1 Institute of Physics and Astronomy, University of Potsdam, 14476 Potsdam 
e-mail: lida@astro. physik. uni-potsdam. de 

2 Massachusetts Institute of Technology, Kavli Institute for Astrophysics and Space Research, 70 Vassar St., Cambridge, MA 02139, 
USA 

3 Leibniz Institute for Astrophysics Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany 

4 Department of Physics and Astronomy, East Tennessee State University, Johnson City, TN 37663, USA 

5 South African Astronomical Observatory, PO Box 9, Observatory 7935, Cape Town, South Africa 

Received September 15, 1996; accepted March 16, 1997 


ABSTRACT 

Context, ji Cep-type variables are early B-type stars that are characterized by oscillations observable in their optical light curves. At 
least one p Cep-variable also shows periodic variability in X-rays. 

Aims. Here we study the X-ray light curves in a sample of p Cep-variables to investigate how common X-ray pulsations are for this 
type of stars. 

Methods. We searched the Chandra and XMM-Newton X-ray archives and selected stars that were observed by these telescopes for 
at least three optical pulsational periods. We retrieved and analyzed the X-ray data for k Sco, p Cru, and a Vir. The X-ray light curves 
of these objects were studied to test for their variability and periodicity. 

Results. While there is a weak indication for X-ray variability in p Cru, we find no statistically significant evidence of X-ray pulsations 
in any of our sample stars. This might be due either to the insufficient data quality or to the physical lack of modulations. New, more 
sensitive observations should settle this question. 
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1. Introduction 


Oscillating stars can be found almost everywhere in the HR 
diagram. The hottest and most massive oscillating stars are 
6 Cephei-type variables. Born with masses between 8 M 0 and 
18 M 0 , and while still young and burning hydrogen in their 
cores, these B0-B2 type stars pulsate with periods of a few 
hours. The physical mechanism that drives these oscillations 
is understood well and attributed to changes in the opacity in¬ 
side the star during the pulsation cycle (“/r-mechanism”, Dziem- 
bowski & Pamiatnykh 1993). 

Like other hot massive stars, p Cep-type variables drive stel¬ 
lar winds by their intense UV radiation. Photons that are scat¬ 
tered or absorbed in spectral lines transfer their momentum and 
thus accelerate the matter to supersonic velocities. This driving 
mechanism is unstable. It is generally believed that the wind 
instability results in shocks where part of the wind material is 
heated to X-ray emitting temperatures (e.g., Feldmeier et al. 
1997a). 

In p Cep-type variables, the deposition of mechanical energy 
from stellar pulsations can provide additional heating of the in¬ 
ner wind regions. This idea was put forward to explain the ob¬ 
served excess in the extreme UV spectrum of p CMa (Cassinelli 


* The scientific results reported in this article are based on obser¬ 
vations made by the Chandra and XMM-Newton X-ray Observatories, 
data obtained from the Chandra and XMM-Newton Data Archives, and 
observations made by the Chandra and XMM-Newton and published 
previously in cited articles 


et al. 1996). Recent work on classical Cepheids (Neilson & 
Lester 2008; Engle et al. 2014) has demonstrated that pulsations 
may power X-ray emission even from such cool stars. Thus, it 
seems reasonable to assume that stellar oscillations are closely 
linked with the X-ray production in p Cep variables. 

Recently, it has been found that X-rays from the strongly 
magnetic p Cep-variable P CMa pulsate in phase with the opti¬ 
cal light curve 1 (Oskinova et al. 2014). To investigate whether 
p CMa is a unique anomaly or if other p Cep-variables are also 
X-ray pulsators, we searched the X-ray archives for observa¬ 
tions of p Cep-variables. The archival data for a selected sample 
of variable were analyzed to extract the X-ray light curves (see 
Figure 1) and study their variability. In this Research Note we 
present the results of this study. 


2. Previous investigations of X-ray pulsations in 
/? Cep-variables 

Despite being a numerous and astrophysically important class of 
objects, p Cep-variables have not been studied well in X-rays. 
The first X-ray survey of six p Cep-type stars was performed 
with the Einstein observatory (Agrawal et al. 1984). No corre¬ 
lation was found between X-ray and pulsational, rotational, or 
binary properties. 

The interest in X-ray properties of p Cep-variables was re¬ 
newed with the launch of the Rosat X-ray observatory. Four 

1 http://sci.esa.int/xmm-newton/541Ql/ 


Article number, page 1 of 5 





A&A proofs: manuscript no. arXbcep 



0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 


Time [hr] 

Fig. 1 . X-ray light curves in 0.2-10.0 keV band of (5 Cep-type variables 
binned by 1000 s. The count rate of f 1 CMa is at least an order of mag¬ 
nitude higher than for the other stars. The horizontal axis shows the time 
after the start of the observation. Vertical bars represent lcr errors. From 
top to bottom: the EPIC PN light curve of p 1 CMa; the HETG (added 
+1 and -1 order) light curve of ft Cru, where the lower red light curve is 
in the hard band (1.0-10.0keV) where the detection of pulsations was 
claimed by Cohen et al. (2008); LETG (added +1 and -1 order) light 
curve of a Vir; EPIC PN light curve of k Sco. 


P Cep stars were among the sample of 12 nearby near-main- 
sequence B-stars observed by Rosat (Cassinelli et al. 1994). Tim¬ 
ing analysis was performed for all targets, but no variability was 
found except those attributed to the spacecraft wobble. 

Rosat observations of p Cep-variables were further analyzed 
for the presence of X-ray pulsations by Cohen (2000), who men¬ 
tions detecting periodic X-ray variability in four /3 Cephei-type 
variables. Specifically, for the star p Cephei, it was shown that 
the hypothesis of a constant source can be rejected with 98% 
confidence (see Fig. 5 in Cohen 2000). Cohen (2000) reports that 
the X-ray variability of [i Cep is periodic with the same period as 
stellar pulsational period in the optical. 

The binary system /I Cep consists of a primary magnetic 
B2III-type star (Donati et al. 2001) and a secondary B6-B8e type 


star (Schnerr et al. 2006). The magnetosphere of the primary was 
studied in detail by Donati et al. (2001). It was predicted that an 
equatorial disk-like structure would be present around this star. 
A signature of such a disk was found using interferometry by 
Nardetto et al. (2011). 

The occultations of X-ray emission sites by the disk should 
lead to observable X-ray variability on the rotational time scale. 
To test this hypothesis, the more sensitive XMM-Newton and 
Chandra observations of p Cep were obtained. However, only 
limited evidence of a modulation in the X-ray emission was 
found (Favata et al. 2009). The variability was observed at a 5% 
level but in anti-phase with the model predictions. Furthermore, 
Favata et al. (2009) performed a dedicated time analysis of the 
XMM-Newton light curves but did not uncover any periodicity, 
thus questioning the results of Cohen (2000). It was concluded 
that “the fact that Cohen (2000) did not perform a period search, 
but rather assumed the period and fitted the amplitude and phase, 
makes it likely that the period reported by Cohen (2000) is spu¬ 
rious” (Favata et al. 2009). 

An assumed period was also fitted to the X-ray light curve 
P Cep-variable /j Cru obtained with the Chandra X-ray observa¬ 
tory (Cohen et al. 2008). These authors applied a Kolmogorov- 
Smirnov test to the unbinned photon arrival times, as well as a 
X~ fitting of a constant source model to the binned light curve 
and found no evidence of X-ray variability. However, when a 
grid of test periods was fitted to the phased photon arrival times, 
it was found that for the X-ray light curve in the hard band 
(hv > 1 keV), the rejection probabilities for some of the assumed 
periods were quite low. This was the basis for suggesting that 
the periodic variability is detected with a period P = 4.588 h, 
which is the primary pulsation mode period in p Cru. However, 
Aerts et al. (2014) give the dominant frequency in p Cru as 
5.95867 d _1 (or 4.0277 h). Polarimetric measurements for p Cru 
did not reveal any magnetic field (Hubrig et al. 2006). We revisit 
the archival X-ray observations of p Cru in Section 3.1. 

Another p Cep-variable that has been studied in X-rays is 
P Cen. This is a binary system consisting of two p Cep-type 
stars of nearly equal mass with an orbital period of 357 days 
(Ausseloos et al. 2006, and references therein). The primary ro¬ 
tates significantly faster (usin; = 190 ± 20kms _1 ) than the sec¬ 
ondary (u sin i = 75 ± 15 km s -1 ), as is evident from its rotation- 
ally broadened lines. The magnetic nature of the secondary, the 
slower rotating component in p Cen is reported by Alecian et al. 
( 2011 ). 

Sensitive observations of p Cen with XMM-Newton are re¬ 
ported by Raassen et al. (2005). The X-ray light curves with 
different time binnings and energy bands were extracted from 
data obtained by all XMM-Newton instruments and their com¬ 
binations. Scargle periodograms were produced but no period¬ 
icities with high statistical confidence found in any of the data 
sets. Moreover, no hint of periodicity with the period seen in op¬ 
tical light curves was found even with low statistical confidence. 
However, when the light curves were compared with the con¬ 
stant source hypothesis, probabilities of 83% (for EPIC PN) and 
49% (for EPIC MOS) were determined for the source not being 
constant. 

Oskinova et al. (2011) considered a small sample of p Cep- 
variables with well known pulsational behavior and existing X- 
ray observations and searched for correlations between pulsa¬ 
tional and X-ray properties. After confirming previous results 
known from the Einstein and Rosat works, no clear relations 
could be found. This could be due to the very limited sample. 

The sample of p Cep-variable with X-ray light curves of suf¬ 
ficient quality is very small. It is observationally challenging to 
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Table 1. XMM-Newton and Chandra observations of p Cep-variables covering at least three stellar pulsational periods 



Sp 

Period A 

[h] 

Obs 

Counts 13 

p-\ 

puls 

< 

X-ray period? 

a d 

max 

Comment^ 

Reference 

P Cen 

B1III 

3.684 

XMM 

25000 

6 

no 

< 3% 

mag 

Raassen et al. (2005) 

P Cep 

B2IIIe 

4.572 

XMM 

2500 

9 

no 

<6% 

mag 

Favata et al. (2009) 

P Cru 

B0.5IV 

4.589 

cxo 

225 f 

4 

no 

< 26 f % 

non-mag 

this work 

k Sco 

B1.5111 

4.795 

XMM 

1700 

3 

no 

< 11% 

? 

this work 

a Vir 

B1.5IV-V 

6.521 

CXO 

470 

3 

no 

<21% 

? 

this work 

CMa 

B0.5-1IV 

5.030 

XMM 

12800 

6 

yes 

10% 

mag 

Oskinova et al. (2014) 


A - main period of stellar oscillations as listed in the “Catalog of Galactic p Cephei stars” (Stankov & Handler 2005); 
B - number of counts per optical pulsation period (EPIC PN for XMM-Newton)\ 

C - number of pulsation periods covered by X-ray light curve; 

D - amplitude of X-ray pulsations or its upper limit for non-detections; 

E - magnetic field presence; 

F - in the hard band where the pulsations were suggested by Cohen et al. (2008); 
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Fig. 2. Power density spectra based on the Fourier transform (Scargle 1982; Horne & Baliunas 1986) of X-ray light curves. The periodogram 
was calculated from v = l/P to v = Ng/P with a spacing P 1 , where Ng is the number of data points. Various false alarm probability levels are 
marked. From left to right: 1) combined zero and first-order Chandra HETGS light curve of pCm in hard (1,-lO.OkeV) band where pulsations 
were suggested by Cohen et al. (2008). There is about a 50% probability that this period is spurious. 2) zero-order Chandra’s HETGS light curve 
of /3Cru in full (0.2-10.OkeV) band, there is more than a 50% probability that this period is spurious; 3) EPIC PN light curve of CMa in hard 
band (0.7-10.0 keV). 


obtain X-ray light curves of high enough quality to establish pe¬ 
riodic variability with high statistical significance. The stellar os¬ 
cillations periods are typically four to six hours. To confidently 
detect pulsations in a light curve, the observations should be long 
enough to cover a number of periods. (In this study we adopt the 
minimum of three periods.) At the same time, the source should 
be X-ray bright in order to have enough counts per time bin for 
a meaningful statistical analysis. In Table 1 we summarize the 
objects that satisfy the above criteria. 

3. Looking for the variability in X-ray light curves of 
[} Cru, k Sco, and a Vir 

3.1. p Cru 

We retrieved the archival Chandra X-ray observations of (i Cru 
(see Table 2) and analyzed them using the latest calibration. The 
X-ray light curves were extracted in different energy bands as 
defined in Cohen et al. (2008): full (0.4-10.0 keV), soft (0.4- 
1.0 keV ), and hard (1.0-10.OkeV). The light curves were binned 
in the 1000 s and 3600 s bins, and examples of light curves are 
shown in Fig. 1. 

For HETG dispersed photon data, we computed light curves 
using the ACIS-Grating-Light-Curve ( aglc ) program (used by 
the Chandra Grating-Data Archive and Catalog processing (TG- 
Cat, Huenemoerder et al. 2011)), in order to also produce a 
curve in the hard band, in addition to the default broad-band 
curve available from TGCat. Zeroth-order curves (not shown) 


were created with the CIAO (Fruscione et al. 2006) program, 
dmextract. The program aglc computes rates with the explicit 
exposure records for each CCD, using the user-specified com¬ 
bination of grating types, orders, and bandpasses. For relatively 
low count-rate systems, such as those analyzed here, with time 
binning comparable to the dither period and with broad energy 
ranges spanning all CCDs, chip-to-chip differences and chip gap 
effects are negligible. Furthermore, for the dispersed photons, 
there is no concern with light leaks, as there could be for zeroth 
orders of optically bright objects. 

Using the^y 2 test, the X-ray light curves were compared with 
those expected from the constant source hypothesis. As an ex¬ 
ample, for the zeroth-order light curve in full band and binned 
by 1000 s, the probability value is P & 0.168, which implies that 
the probability that the observed variability is solely a chance oc¬ 
currence is less than 17%, which is statistically insignificant. 
Thus, the null-hypothesis assuming that the source is constant 
cannot be rejected. On the other hand, for the same light curve 
but binned into 3000 s, P ~ 0.0044, which is statistically signif¬ 
icant at the level a = 0.01. In the hard band, for the combined 
zeroth- and first-order light curve, P ~ 0.06 is statistically in¬ 
significant with a = 0.05. 

Thus, similar to the case of ft Cen, there are marginal indica¬ 
tions of source variability. To investigate the presence of periodic 
modulations further, we constructed Scargle periodograms, with 
some examples shown in Fig. 2. A peak with a false alarm proba¬ 
bility of about 50% is present at the frequency close to the stellar 
pulsational period. However, this is not the most prominent fea- 
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Table 2. Log of analysed observations 


Object 

Telescope 

ObsID 

Start time 

Duration 

[ks] 

Count rate 
[s- 1 ] 

/SCru 

CXO/HETG 

2575 

2002-02-01 

75 

0.03 (-1 + 1 order) 

k Sco 

XMM 

0503500201 

2008-03-12 

64 

0.1 (pn) 

a Vir 

CXO/LETG 

4509 

2004-03-09 

90 

0.02 (-1 + 1 order) 


ture in the periodograms. For comparison, we show in Fig. 2 a 
periodogram of the £ ;l CMa in 0.7-10.0keV (hard) band, where 
the peak at the stellar pulsation frequency has a false alarm prob¬ 
ability less than 1%. 

Thus, from reanalysis of the Chandra data on /3 Cru, we must 
conclude that in the statistical sense, it is likely that the period 
reported by Cohen et al. (2008) is spurious. Better quality obser¬ 
vations are needed to establish whether the periodicity is indeed 
present in /3 Cru’s X-ray light curve. 

3.2. X-ray variability of k Sco 

The spectroscopic binary system, k Sco, has two B-type com¬ 
ponents orbiting with a period of 195 d. The primary (B1.5111) is 
responsible for the pulsational variability (Harmanec et al. 2004) 
observed photometrically with a main period of 4.797 h (Handler 
& Schwarzenberg-Czerny 2013). 

We retrieved the archival XMM-Newton observations of k 
Sco (see Table 2) and reprocessed them using the most recent 
data analysis software (SAS 14.0.0) and calibration files. The 
source k Sco has a visual magnitude my = 2.4, therefore the 
thickest optical light blocking filter was used in observation. For 
the XMM-Newton detectors operating in full-frame mode and 
with the thick filter, no optical loading is expected for stars with 
my > -2 mag (for the PN camera). 

Unfortunately, the observations were affected by the high 
background from proton flares. The observations were filtered 
using good time intervals resulting in a useful exposure time of 
44 ks. Here, k Sco is sufficiently isolated without nearby X-ray 
sources. From the region with radius 25" and centered on the 
coordinates of k Sco, we extracted the EPIC PN light curve in 
0.2-10.0keV energy band and restricting the patterns to single 
and double events. Background regions were carefully selected 
on the same CCD. The net light curves with different time bin- 
nings were produced using the epiclccorr task. As a next step, 
the X-ray light curves were tested for variability. As an example, 
in the full band and binned with a 1000 s light curve, the proba¬ 
bility value P » 0.83 is not significant even with a = 0.1. Thus, 
the null hypothesis about the X-ray light curve from k Sco be¬ 
ing constant cannot be rejected. No useful information about the 
presence of periodic variability could be obtained on the basis of 
Fourier periodograms, owing to the limited quality of the data. 

3.3. X-ray variability of a Vir 

Spica (a Vir) is a spectroscopic binary with a period of ~ 4 d in 
an eccentric orbit. The primary component (B1.5IV-V) is classi¬ 
fied as a /? Cep-type variable, but its pulsational behavior is not 
stable (see discussion in Palate et al. 2013). 

Spica was observed by Chandra LETG for 90180s in 2005- 
03-15 ( Chandra observation identifier 4509). We retrieved the 
X-ray light curves (see Fig. 1) from the TGCat and tested them 
for variability. As mentioned above, there is no concern about 



1 IP [ 10‘ 4 s' 1 ] 

Fig. 3. Power density spectrum based on the Fourier transform (Scargle 
1982; Horne & Baliunas 1986) of X-ray light curve (0.4-10.0 keV) of a 
Vir. The false alarm probability at 50% level is marked. 

optical or UV light leakage since we are not using the zeroth- 
order photons. For extremely high UV flux sources, the LETG 
coarse support structure acts as a crude disperser, and it can 
sometimes produce spurious signal for LETG used with HRC-S 
(Drake 2010), but we see no such evidence in this LETG/ACIS-S 
spectrum. 

Based on the ;y 2 -test, the chance probability value P ~ 0.19 
is not significant. We went on to conduct a period search, but all 
peaks in the Fourier periodogram are statistically insignificant 
(see Fig. 3). No useful information on the X-ray variability was 
obtained from other statistical tests either. Moreover, in the case 
of a Vir with its unstable pulsational behavior in order to reach 
any conclusions about the correlation between stellar oscillations 
and X-ray variability, it would be necessary to observe the star 
simultaneously in X-rays and in optical. 

The Chandra observations of a Vir lasted for about 25 h, thus 
covering about a quarter of the orbital period of this binary star. 
Although no optical eclipses have been observed in this system 
(Palate et al. 2013), an influence of binarity on the X-ray light 
curve cannot be excluded a priori. 

4. Discussion 

Our study of a small sample of X-ray lightmcurves from fi Cep- 
variables did not reveal any periodic variability. However, the 
X-ray count rate of the only star where the X-ray pulsations are 
clearly detected, f 1 CMa, is much greater than for the other stars 
considered here. It is important to assess what the upper limit is 
on an amplitude of modulation that could have been detected in 
our sample stars with the available data. 

We therefore wish to study a star with a periodic signal of 
relative amplitude A and N counts per pulsational period. The 
noise level can be approximated as 2 y/N/n p , where n p is the 
number of periods. We assume that a significant detection must 
be a factor 4x higher than the noise level. Thus, we can derive 
the minimum detectable amplitude for a given observation. The 
results are shown in Table 1. Only for k Sco the sensitivity of 
observations was marginally sufficient to detect a variability on 
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a level of about 10%. The data sets for [i Cru and a Vir are not 
good enough to uncover periodic variability on a level similar to 
that seen in CMa even if it were present. 

While we did not found any periodic pulsations, our analysis 
showed that some of the sample stars show marginal evidence of 
X-ray variability. For some /j Cep-type variables, such variability 
has also been previously reported in the literature. 

Significant stochastic X-ray variability is expected in the 
models of X-ray generation in stellar winds (Feldmeier et al. 
1997b; Oskinova et al. 2001b). However, the most sensitive X- 
ray observations of an O-type star (the XMM-Newton observa¬ 
tions of £ Puppis with the total exposure time of about 1 Ms) re¬ 
vealed a surprising lack of short time-scale stochastic variability 
(Naze et al. 2013). 

On the other hand, the evidence is accumulating that the X- 
ray variability associated with stellar rotation may be present in 
massive star winds. Berghoefer et al. (1996) observed the super¬ 
giant £ Puppis for 11 days and found modulations with a fre¬ 
quency of 1.44 d 1 in the Ha line profiles, as well as in the X- 
rays in the 0.9-2.0keV band. They reported an amplitude of X- 
ray variability of 6% and suggested that these variations are pe¬ 
riodic. This suggestion was not confirmed by later studies based 
on ASCA and XMM-Newton observations (Oskinova et al. 2001a; 
Naze et al. 2013). Howarth & Stevens (2014) reported the detec¬ 
tion of a period P = 1.780d in optical light curves of £ Pup¬ 
pis and suggested that this period may be due stellar pulsations. 
While indeed the X-ray modulations on a time scale of a day are 
present in £ Puppis, no coherent periodicity in X-ray emission 
has been confirmed for this object so far. 

Periodic modulations of X-ray emission were reported for 
the O-type dwarf £ Oph (Oskinova et al. 2001a). However, also 
in this case, we lack any information on the repeatability of X- 
ray modulations. Two other O-type stars, £ Per (Massa et al. 
2014) and £ Ori (A. Pollock, private communication), also show 
modulations of X-ray emission on a time scale of days. In¬ 
terestingly, XMM-Newton observations of a Wolf-Rayet type 
star, WR6 showed that its X-ray light curve changes quasi- 
periodically (Ignace et al. 2013). We speculate that, similar to 
WR6, in the OB-type stars the X-ray variability owing to the 
rotation has quasi-periodic character. 

A small fraction of B-type stars possess strong magnetic 
fields (e.g., Hubrig et al. 2006; Morel et al. 2014, and refer¬ 
ences therein). If a magnetic field has a dipole configuration, it 
can channel the wind toward the magnetic equator, where wind 
streams from opposite hemispheres collide and produce a strong 
shock (magnetically confined wind shock (MCWS), Babel & 
Montmerle 1997; ud-Doula et al. 2014). The MCWS predicts 
the X-ray variability of magnetic massive stars associated with 
the stellar rotation due to the occupation of X-ray emitting re¬ 
gions by the cold torus of matter accumulated in the equato¬ 
rial plane (Donati et al. 2001, 2002). Such X-ray modulations 
on a rotational time scale are sometimes observed (e.g., Gagne 
et al. 1997; Pillitteri et al. 2014); however, this is not ubiquitous 
behavior among magnetic massive stars - some of these well- 
known objects, such as r Sco, do not show any evidence of X-ray 
variability over the stellar rotation period (Ignace et al. 2010). 

The data sets that we consider in this paper do not cover long 
enough time intervals to test the X-ray light curves for modula¬ 
tions due to stellar rotation or orbital motion. Thus, it is not clear 
what the true nature is of X-ray variability observed in some of 
our sample stars. 

However, in none of our sample stars, [5 Cru, k Sco, or a 
Vir, do we find statistically significant evidence for periodic X- 
ray variability. Up to now, the unambiguous presence of X-ray 


pulsations with the same period and in phase with the optical 
light curve has only been detected in the magnetic star £ ] CMa. 
It remains to be seen whether this object is exceptional or if new, 
better quality observations will reveal similar properties in other 
/? Cep-variables. 
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